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1 Introduction 



All quantum mechanical models of the infinite sequences A n , B n ,C n , BC n , D n 
and in addition some exceptional sequence models which were defined and 
proved to be completely integrable by Olshanetsky and Perelomov by the 
"Hamiltonian reduction method" [0J possess limiting versions of Sutherland 
and Calogero type that can be shown to be exactly solvable by a simple 
trick. Namely, their Schrodinger operators can be transcribed into quadratic 
polynomials of related "hidden" Lie algebras. The spectrum and the eigen- 
functions can then be determined from representation theory, e.g. using 
polynomial spaces . 

This program was formulated in || and successfully applied first to the 
A n sequence in |3j . Then it was carried over to the other sequences and even 
the super symmetric generalizations in |fj. |5[. 

Our aim is to turn the arguments around and to develop an algorithm 
which leads to new exactly solvable models. First investigations were pre- 
sented in || (hereafter quoted as (I)). The program contains two major and 
separate issues: to render a second order differential operator curvature free 
and to find a first order differential operator satisfying an integrability con- 
straint. If both constraints are fulfilled, the whole operator can be trans- 
formed into a Schrodinger operator with a standard Laplacian as kinetic 
energy and a real potential as potential energy terms. 

In this work we apply this method to four-particle models. We start 
from the second order differential operator for A% known from || (Calogero 
case) and || (Sutherland case). After a change of the variable t% to A3 = r| 
we test out all possible first order differential operators. This way we find 
a new solution leading to a four-particle potential which is automatically 
translation-invariant. In the Calogero case this model was known before J/J, 
but the Sutherland case seems to be new. 



2 The program 

We are interested here in the bound state spectrum of Schrodinger operators. 
The whole analysis is therefore performed in real spaces. Consider a flag of 
polynomial spaces Vn, N e Z> 

V N = span{xi 1 X2 2 ...x^ n |ripi + r 2 p 2 + ••• + r n p n < N} (1) 

(Pi e IN) 



We consider differential operators of first order 

dx a 



D\L = ^4r (2) 



2 



(a a multi-exponent) 
and of second order 



n(2) _ JS] d2 {Vi 
L >[S;a,b]- X dXa Q Xb ^ 

that leave each space Vn invariant. If 

p=(l,l,...,l) (4) 

then the operators (0) generate the full linear (inhomogeneous) group of H n 
and the operators of second order (|3]) can be obtained as products from the 
first order operators, i.e. in (|2|) 

a = e<*>, et ] =5 a h (5) 

and in @ 

5 = e (a) + e (6) ( 6 ) 

Now we consider a candidate for a future Schrodinger operator 

D = 9{S;a,b] D {S;a,b] 
ct,a,b 

The eigenvectors and values of D in V^r can be calculated easily by finite 
linear algebra methods (the number of eigenvector in Vn may be smaller than 
dim V^, see (I)). If we want completely integrable models we must make sure 
that a complete set of involutive differential operators exists. For this task 
Lie algebraic methods may be very helpful. 

The first step in transforming D into a Schrodinger operator is to write 
it symmetrically 

fl = -E^w|- + EU^ (8) 

where 

9ab=J29[a;aM^ ] (9) 

a 

We write because this is the inverse of a Riemann tensor. The Riemann 
tensor g ab is assumed to be curvature free. The task to make it so will not 
arise in this work. But we mention that we developed a minimal algorithm 
to solve this issue. 

Following the notations of (I) we "gauge" the polynomial eigenfunctions 
of D by 

$(x) = e- x{x) ip(x) (10) 
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so that 

e-^ + * = -^E^(v«)^ + ^) (ID 

(g = (det 9 -')-i). 

This is possible if and only if 

Ea*(*)^[ 2 x-kvfl = r.0»0 (12) 

which implies integrability constraints on the functions {r a (x)}. If they are 
fulfilled we obtain a "prepotential" 

p = lnP (13) 

so that 

p = 2 X -\n^ (14) 
In all cases studied, we found solutions for p as follows. Let 



detg- 1 {x) = HP i {x) (15) 

i=l 

where {Pi(x)} are different real polynomials. Then 

r 

P(z) =E7ilnPi(x) (16) 



i=l 



with free parameters ji solves the requirement that {r a (x)} (|12"D belong to 
differential operators leaving each Vjv invariant. In particular 



1 _ _ lf .dP, 



Pi(x) V " ^ 



are polynomials, and each power [a] is bounded by the requirement that the 
spaces Vtv(1) are invariant. Inserting fll5|), ([16D in ([14]) we obtain finally 

X = ^E(7i-i)^(a?) (18) 

i=l 

We will later see that in the case of the models of Calogero type a term 

7olnP (^) (19) 

can be added to p, where 

P {x) = e* 1 (20) 

is not contained in det g~ l as a factor. This prepotential gives rise to the 
oscillator potential. 
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Finally we mention that e~ x is the ground state wave function of the 
Schrodinger operator, as follows from ([Top. 

The expression (I, 6.17) for the potential W(x) contains a term linear in 

Each divergence 

E^)(x) = C« (22) 

ought to be a constant to render the operator (^) formally self-adjoint after 
subtraction of a constant. From now on we shall dismiss all constant terms 
in W(x). 

We can then write the potential as 

in.r) Y^'.jlhj (23) 

= (24) 

lv = \(^-\) (U/O). (25) 
In the cases of this article 

Rij = const if i ^ j (26) 



If we then set 



we obtain 



H = -Vi + \ (*^0) (27) 



with 



W(aO = $>«fl« (28) 

i=l 

T« = - 1) (29) 

3 A model of Calogero type 

In the Schrodinger operator for the A3 Calogero model 

#cal = \ £ + + 51 ^(^ - ^)" 2 (30) 



5 



we separate the centre-of-mass part by introducing relative coordiantes 

X = X\ + %2 + x 3 + £4 
Vi = Xi-^X, ie {1,2,3} (31) 

Following Q we introduce elementary symmetric polynomials as new coor- 
dinates instead of the {r^} 

l<i<j<4 

T 3 = ViVjVk 
l<i<j<fe<4 

ta = ymyzVA (32) 

Using the ground state wave function 

%(x) = ^^exp j-^^j (33) 

with V\ any root of 

g 1 = v x {v 1 -\) (34) 
and V(x) the Vandermonde determinant, we gauge H ca \ to 

#cal = * (^)~ 1 ^cal*o(2:) 

= #c.o.m. + H rc i + const (35) 
Then (see || and I (7.1) - (7.6) with s = 0, W2 = u> 3 = u> 4 = 1) we get 

2# rcl - - ^ <u 

+ first order differential operator (36) 

with 

/ -2r 2 -3r 3 -4r 4 \ 

<T 1 = -3r 3 -4T 4 + r| + ±t 2 t 3 (37) 
V -4r 4 + ±t 2 t 3 -2r 2 r 4 + |r 3 2 / 

The first order differential operator is discarded. 

The relative motion Hamiltonian ( |36"D reflects the symmetry of H ca \ ( |30| ) 
under 

.Xi -> -Xi (38) 

namely it is symmetric under 

r 2 -> r 2 , r 3 -> -r 3 , r 4 -> r 4 (39) 
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Thus it leaves polynomial spaces that are even or odd in r 3 separately in- 
variant. Now we reduce these spaces and consider polynomials only that are 
even in r 3 : 

r 2 = A 2 , r 2 = A 3 , r 4 = A 4 (40) 

9 -24 9 



dr 3 6 d\ 3 

d 2 d 2 d 

sf = ^m +2 ^ (41) 

(this substitution fl4"0| ) was invented in j|, Q). Multiplying the second row 
and column with 2A 3 we obtain 

9+ -» Gab (42) 
Whereas the quadratic differential operator (K6l) with g~ h as in d37|) and 



P^t) = -4det0 _1 (I, (7.22), and appendix (A.l)) (43) 

Po(r) = e T2 (44) 

r (1) = (-12,0,-2r 2 ) (45) 

r (0) = (-2r 2 ,-3r 3 ,-4r 4 ) (46) 

allows to reconstruct H c &\ completely, we have now an additional factor in 

detG- 1 = 4A 3 (-^)Pi(A) (47) 

(Pi (A) is obtained from P\{t) (A.l) by substitution (|40|)). 
The polynomial spaces ([!]) are now such that 

n = 3, px = p 3 = 1, p 2 = 2 (48) 

For the first order differential operators (^) we have 

[a; a] G j [0; a] all a; 

[ e W;l], [ e W;3], og{1,3}; 
[e (a) ;2], alia; 

[e (a) + e (fe) ;2], a, G {1,3}| (49) 

If we multiply two of these we obtain possible second order differential oper- 
ators (0), but in addition we have 

0;a,b] G {[e {2) ;a,b] a, b G {1,3}} (50) 
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Now the program of the preceding section is set in action: From 

Pi (A) = -^detG- 1 (51) 
A3 

P 2 (A) = A 3 (52) 

Po(A) = e M (53) 



we obtain 



r (1) = (-12,0,-2A 2 ) (54) 
= (- 6 ,4(A^-4A 4 ),A 2 ) (55) 
r<°> = (-2A 2 ,-6A 3 ,-4A 4 ) (56) 



and 



16 

R n = — [2X1 + 16A^A 4 - 96A 2 A^ + 9A^A 3 + 108A 3 A 4 ] (57) 
Pi 

P22 = ^[A^-4A 4 ] (58) 
^3 

Poo = -2A 2 (59) 

Pn leads us back to the Calogero and Poo to the oscillator potential. The 
new potential following from P 22 is obtained as follows. We factorize A3 
(eliminating y A = -y x - y 2 - 2/3) 

A3 = A = ( Vl + y 2 f{ Vl + y 3 ) 2 (y 2 + y 3 ) 2 (60) 

and perform a fractional decomposition 

Y= E (Vi + Vj)-' (61) 

3 l<«<i<3 



We find 
so that 



Q = \ 2 2 - 4A 4 (62) 
P22 = 16 ( x * + x i - x k-xi)~ 2 (63) 



3 inclcp. terms 

Using (|29|) we obtain finally the new model 

^ = ^g(-^ 2 +^ 2 



+ 9i E 

l<i<j<4 

+ 92 (xi + Xj -x k - xi)~ 2 (64) 



3 indcp. terms 



with 



g 2 = 2u 2 (u 2 - 1) 



(65) 



In the course of this derivation we have reduced the polynomial spaces of 
the standard A% Calogero model to spaces of polynomials that are even in 
T3. But afterwards we multiply them with 

*t = T? (r 3 >0) (66) 

So for v 2 = we have recovered the even and for v 2 — 1 the odd spaces. In 
these cases g 2 = as shown by d65|). Thus the new potential arises by an 
interpolation between even and odd spaces. 



4 A model of Sutherland type 

Having constructed a model of Calogero type the existence of a corresponding 
model of Sutherland type can be expected. We start from our formulas for 
the N = 4 Sutherland model in (I, Section 7) and put everywhere 

Wi = l, i G {2, 3, 4}, s = 2 (67) 



Again we change variables as in (f40|) , ( f41|) and obtain 



det(G- 1 )=4A 3 (-i)P 1 (A) (68) 



as in (|47|) , and 

P 2 (A) = A 3 (69) 

as in (|52|). Pi (A) is given as in the appendix (A. 2). There is no place for a 
P (A). The reader is advised to try out what happens with a -Po(^)- 
Calculations become very lengthy. We find 

r« = (2(a 7 -l)(a 7 -3)A 4 -2(a 7 + 5)A 2 -12, 
-24A 3 , 

+2(a 7 - 9)A 4 - 2A 2 ) (70) 



rf = (-(a 7 -l)(a 7 -3)A 4 + (a 7 -7)A 2 -8, 
+4[(a 7 - 1)A 4 - A 2 ] 2 - 8A 3 - I6A4, 
-(a 7 + 3)A 4 + A 2 ) (71) 

The parameter a 7 appears only in the Lie algebraic formulation of the same 
(a 7 independent) Sutherland model. 
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The potentials Rij come out as 

Rn = ^-{Ri} (72) 
with Ri given in the appendix (A. 3) 

R22 = ^{(a 7 -l)A 4 -A 2 ] 2 -4A 4 } 

+ irrelevant terms (73) 

and R12 is constant (irrelevant). We introduce (see (I), eqs. (7.16)-(7.18)) 
the shorthand for the trigonometric ^-symmetric functions 

a 2 = sin Vi sin Vj cos Vk cos yi 



-- sin 2 (^ + % ) (74) 

i<«<i<3 



a 3 = sin Vi sin Vi sin Vk cos Vi 

= - ]J sin (yi + Vj) ( 75 ) 

i<«<i<3 

<j 4 = sin y 1 sin ?/ 2 sin y 3 sin y 4 (76) 



Since 



A 3 = a\ (77) 
we study the fractional decomposition 

Y= £ (sm^ + %))- 2 (78) 

3 1<*<?'<3 

For Q we find 

Q = (a 2 - 2a 4 ) 2 - 4er 4 + <x 2 (79) 

On the other hand we get by insertion of A 2 , A 4 (£ = 2A 2 , ?? 4 = 4A 4 in I (7.9) 
and (7.11)) 

R22 = ^ + const (80) 
A3 

Thus we end up with a Sutherland model 

H suth = \Y{-^2+^) 

+gi Y ( sin ( x i ~ x j))~ 2 

l<«<i<4 

+ j92 Y ( Sil1 \( X i + X 3 ~ X k ~ X l)Y 2 ( 81 ) 

3 cases 



with g 2 as in (|65|). The remarks made at the end of the preceding section 



remain true in this case. 
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5 Concluding remarks 



Our technique can also be applied to the A 2 model whereas the cases Ajv-i, 
N > 4, cause hitherto unsolved problems. In the latter cases the form of the 
trigonometric models in the variables {u n } 2 analogous to ( |74|) - (|76|) is also 
known ||. The trigonometric variables {^n}^ -1 introduced in || in which 
all A N _i models were expressed, are algebraically related with the {(J n } 2 but 
the explicit form of this relation is not known. The trigonometric variables 
introduced in || eqns. (2.17), (2.18) are $lr)i,$sr)i respectively, those of 0, 
eqns. (3.4), (3.5) are still different. 

The rational A 2 and G 2 models are equivalent if 

u(A 2 ) = u(G 2 ) + fi(G 2 ) (82) 

where the respective coupling constants are (see 0) 

g(A 2 ) = u(A 2 )(u(A 2 ) - 1) (83) 
g(G 2 ) = v{G 2 ){v{G 2 )-l) (84) 
9l (G 2 ) = 3MG 2 )(MG9-1) (85) 

Indeed, their Lie algebraic version (0), eqns. (2.8), (4.5)) are identical if 
([82] ) is satisfied. The trigonometric versions are, however, inequivalent. 

Applying our technique to the A 2 model gives the G 2 model, both in the 
rational and the trigonometric cases. Since the A 2 model is obtained from the 
G 2 model by specialization of one coupling constant, we should denote both 
models together as AG 2 . The extension of the A3 model found by us could to 
be denoted AG% correspondingly. The Lie algebras A 3 and D 3 are identical. 
But A 3 and D4 models involve four particles, the center-of-momentum motion 
of the -D 4 model is not separable. Its potential is 

g \ v ( x i ~ x j) + v ( x i + x i)\ ( 86 ) 

l<i<j<4 

whereas our AGs model has potential 

\giv(xi - xj) + jg 2 v(vi + yA) (87) 

They are obviously related, but ( [37|) differs from (|86D by independence of the 
coupling constants and translational invariance. 

Finally we want to mention that we have derived the AG% model also by 
the use of the coordinates /ii, fi 2 , ^ 3 from 

/"1 = Y.y} = -lT 2 (88) 

1=1 
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^4 



E2 2 
ViVj 

l<i<j<4 

E2 2 2 
ViVjVk 

l<i<j<k<4 
2 2 2 2 2 



+2r 4 + r 2 2 



-2r 2 r 4 + To 



so that 



^4 = ^2 



1 



(89) 
(90) 
(91) 

(92) 



This derivation goes along the same lines but is more laborious. It is re- 
markable that in this case one has to use an infinite ensemble of flags of 
polynomial spaces 



V N = span-JX 1 /^ 2 /^ 3 ; raipi + n 2 p 2 + n 3 p 3 < N} 



(93) 



namely one flag for each admissible triplet (pi,P2,Ps) £ IN 3 - In the following 
table we give some of these triplets. The minimal admissible triplet is (2, 3, 5). 



(1,2,4) 


(2,3,5) 


(3,5,10) 


(1,3,6) 


(2,3,6) 


(3,7,14) 


(1,4,8) 


(2,4,5) 


(3,8,16) 


(1,5,10) 


(2,5,6) 


(3,10,20) 


(1,6,12) 


(2,5,10) 


(3,11,22) 




(2,6,7) 




(l,P2,2p 2 ) 




(3,P2,2p 2 ) 



Table 1: Table of admissible triplets (pi,p 2 ,P3)- 

Criterion for admissibility is that the differential operator D (7), (8) maps V N 
in Vn and possesses nontrivial eigenvectors. 



6 Note added 

Recently we were able to prove that the translation invariant four-particle 
model AG 3 considered here is equivalent to the translation non-invariant 
three-particle model belong to B 3 . The arguments rely on a discussion of 
the underlying Coxeter groups and their invariants. 
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7 Appendix: Some formulas 

In the case of the Calogero A 3 model we have from (1,(7.22)) 

P^t) = 27r 3 4 - 2567-f + 128r 2 2 r 4 2 - 16r 2 4 r 4 + 4r 2 3 r 3 2 - 144r 2 r 3 2 r4 (A.l) 

In the case of the Sutherland A 3 model we have (we give these functions for 
a 7 = 3 only) 

P X (A) = -64A 3 A 3 A 4 - 16A^A 4 - 48A 2 A 4 A 2 - 144A 2 A 4 A 3 + 528A 3 A 2 A 2 + 256A 3 A 2 A 2 
-384A 3 A 2 A^ + 128A^A^ - 768A 2 A^ + 27\l + 192A 3 A^ - 48A 4 A 3 ! + 8A^ 
+48A 3 ! A 4 ; - 384A 3 A^ + 192A 3 A 4 - 256A 3 + 768A^ + 4A^ - 208A 2 : A 3 A 4 (A.2) 
+36A 2 A 2 - 768A^ + 256A^ - 16A 4 A 4 + 256A 2 i A 2 - 1024A 2 A 4 ; + 144A 4 A 2 
-512A^ + 896X 2 2 \t - 768A 2 A^ + 1536A 2 A 4 + 4A 3 A^ + 4A 3 A^ 

and 

Ri = 272A 3 A 2 i A 4 + 108A 4 A 3 + 56A 2 'A 4 + 216A 2 A 4 A^ + 9XjX 3 + 732A 2 A 4 A 3 + 2X^ 
-2352A 3 A 2 A^ - 1136A 3 A^A^ + 1728A 3 A 2 A^ - 96A 2 A^ + 16A^A 4 - 736A 2 ,A 4 ! 
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+3648A 2 A^ - 108A^ - 972A 3 A^ + 216A 4 A^ - 34A^A 3 ! - 216X 2 3 \l (A.3) 
+1728A 3 A 4 ; - 864A 3 A 4 ; + 1152A| - 3456A 4 L - 18A^ + 2X 6 2 + 896A^A 3 A 4 
-150A 2 A^ + 3456A^ - 1152A° + 72A 2 t A 4 - 1200A^ + 4736A^A^ 
-600A^ + 2240A^ - 4000A^A^ + 3456A 2 A^ - 7008A 2 A 4 = - 4A 3 A^ 
-14A 3 A* 
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